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Computational/Experimental Aeroheating 
Predictions for X-33 Phase II Vehicle 


H. Harris Hamilton II, K. James Woilmm'nstrr. 
Thomas J, Horvath, and Scott A. Berry t 
NASA Langley Research Center, Hampton. VA 22681 


Laminar and turbulent heating-rate calculations from an “engineering” code and 
laminar calculations from a “benchmark” Navier-Stokes code are compared with experi- 
mental wind-tunnel data obtained on several candidate configurations for the X-33 Phase 
II flight vehicle. The experimental data were obtained at a Mach number of 6 and a 
freestream Reynolds number ranging from 1 to 8 x 10 h /ft. Comparisons are presented 
along the windward symmetry plane and in a circumferential direction around the body 
at several axial stations at angles of attack from 20 to 40 deg. The experimental results 
include both laminar and turbulent flow. For the highest angle of attack some of the 
measured heating data exhibited a “non-laminar” behavior which caused the heating to 
increase above the laminar level long before “classical” transition to turbulent flow was 
observed. This trend was not observed at the lower angles of attack. When the flow 
was laminar, both codes predicted the heating along the windward symmetry plane rea- 
sonably well but under-predicted the heating in the chine region. When the flow was 
turbulent the LATCH code accurately predicted the measured heating rates. Both codes 
were used to calculate heating rates over the X-33 vehicle at the peak heating point on 
the design trajectory and they were found to be in very good agreement over most of the 
vehicle windward surface. 


Nomenclature 

Alt altitude, ft 

h heat transfer coefficient., BTU — sec/ ft 4 

H enthalpy, ft 2 /see 2 

L vehicle length, inches 

M Mach number 

p pressure, lb/ ft 2 

R( i Reynolds number based on length 
R n nose radius, inches 

q heat transfer rate, BTU / ft 2 — sec 

T temperature, (IcyF 

\ velocity, ft /sec 

x. y , z Cartesian body coordinates, inches 

(\ angle of attack, deg 

f surface emissivity 

p density, ships/ ft 2 

Subscripts 

FR Fay and Riddell 

w Wall 

dc Freestream 

Introduction 

T HE Access to Space Study 1 by NASA recom- 
mended the development of a heavy-lift, fully 
reusable launch vehicle (RLV) 2,:J to provide a next- 
generation reusable launch system capable of reliably 
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serving national space transportation needs at greatly 
reduced cost. As part, of this program, the X-33 serves 
as a technology demonstrator for the R.LY. It is envi- 
sioned to be approximately a one-half-scale prototype' 
of the RLY that will demonstrate the technologies re- 
quired to develop and operate 1 a full scale vehicle. 

At the completion of a X-33 Phase I competition. 
Lockheed Martin Skunk Works was awarded the Phase 1 
II contract for the design, development and construc- 
tion of the X-33 flight vehicle. In addition. Lockheed 
is in charge of the flight te'st program. It is an in- 
dustry led effort in partnership with NASA to support 
the design and dewedopmemt work on the 1 X-33 through 
formal task agreements. 

In the 1 design of any hypersonic vediicle. aero- 
dynamic heating is an important issue 1 . This paper 
presents results generated in support of task agree- 
ments with Lockheed to provide experimental aerody- 
namic heating data anel computational fluid dynamic 
(CFD) calculations in support of X-33 aorothormody- 
nainic development, and design. 4 It should be 1 viewed 
as a report of work-in-progress since additional weak 
will be performed during the next several months that 
will also support the final design. 

Results from experimental heating measurements 
and CFD computations are present eel and compared. 
In addition, ongoing efforts in support of these task 
agreements will be discussed. Some of the 1 e*xperi- 
mental heating results presented in this paper were 
in support of the task to study laminar-t.o-turbulent 
boundary-layer transition/' 
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X-33 Geometry 

The X-33 configuration is a lifting body design 
with symmetrical canted fins, twin vertical tails, and 
two body flaps located at the rear of the fuselage. The 
vehicle is powered by a linear Aero-spike engine. It has 
a length of approximately 75G inches from the nose to 
the end of the engine block and a span of approxi- 
mately 920 inches from wing tip to wing tip. The 
geometry of the X-33 vehicle has evolved during the 
course* of both the Phase I and II efforts, 4 thus results 
from several different vehicle geometries (D Loft, Rev 
C, and Rev F) are discussed in this paper. The D Loft 
was the configuration that emerged from Phase I and 
was used in the early part of Phase II. 

The Rev C configuration is a revision of the D- 
loft. For Rev-C, the nose region had a slightly different 
shape to simplify the construction of the metallic pan- 
els used in the Thermal Protection System (TPS) arid 
there were also changes in the base region. Both a top 
and bottom perspective view of this configuration are 
shown in Figs. 1 and 2. respectively. 



Fig. 1 Top view of X-33 Rev— C configuration 



Fig. 2 Bottom view of X-33 Rev— C configuration 

The dihedral of the canted-fin oil both the D -Loft and 
Rev C configurations was 37 degrees. 

The Rev F is a further revision of the vehicle 


geometry to improve 1 its aerodynamic characteristics. 
This configuration has the saint* forebody shape as 
Rev C but the dihedral of the canted fin was low- 
ered to 20 degrees to improve pitch-trim characteristics 
across the speed range and the size of the leeside ver- 
tical tails was increased to improve lateral-directional 
stability at low speeds. Although then 1 have* been ad- 
ditional vehicle modifications since Rev F, these mod- 
ifications have been relatively minor. 

Experimental Method and Test 
Conditions 

Experimental results presented in this paper 
were obtained in the NASA Langley Research Center 
(LaRC) 20-Inch Mach G Tunnel. The tunnel is a hy- 
personic: blow down facility which uses heated, dried, 
filtered air as the test medium. Typical operating 
conditions for this tunnel cover a range of stagnation 
pressures from 30 to 500 psi, stagnation temperatures 
from 760° R to 1000° R, arid freestrearn unit Reynolds 
numbers from approximately 0.5 to 8 x 10 ( * per foot. 
For these operating conditions, the freestrearn Mach 
number ranges between approximately 5.8 and 6.1. 
The tunnel has a two-dimensional nozzle and a test 
section that is 20.5 by 20 inches. A more detailed de- 
scription of the facility is given by Micol. 0 

Experimental heating tests were conducted oil a 
10-inch forebody of the D Loft configuration, repre- 
senting approximately the forward 70 percent of the 
full vehicle, and on a Rev C, 10-inch “tip-to-tail” con- 
figuration in the Langley 20-Inch Mach G Tunnel. The 
tests on the D Loft forebody model were part of a 
study designed to obtain “smooth body or natural” 
boundary-layer transition on the X-33 vehicle. 5 Be- 
cause of the fast paced nature of the X-33 program, 
the models used in the experimental tests presented in 
this paper were not measured to determine how well 
they agreed with the vehicle shape. However, for all 
future experimental heating models, quality assurance 
measurements are planned as a part of the model con- 
struction process. 

The experimental test conditions for these exper- 
iments are given in Tables 1 and 2, respectively. Two 
identical models (designated A and B) were used dur- 
ing the D Loft forebody tests to look at the effect of 
model construction on both the heat transfer and tran- 
sition results. 

All the heat transfer measurements were per- 
formed using a two-color, relative-intensity phosphor 
thermography technique to enable optical acquisition 
of the test data.' 8 With this technique, silica ce- 
ramic wind-tunnel models are slip casted 9 and coated 
with a mixture of phosphors which fluoresce in two 
regions of the visible spectrum (red and green) when 
illuminated with ultraviolet light. The fluorescence in- 
tensity is dependent on the amount of incident light 
and local temperature of the phosphors. By acquiring 
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Table 1 Test conditions for D Loft forebody. 


Run 

Model 

M-* 

n 

(log 

I>3C 

lb/ft- 

Toe. 

dog R 

R <7. 
x 10~“ 

G 

A 

5.95 

40 

12.22 

112.9 

1.8G 

8 

A 

5.97 

40 

23.77 

111.8 

3.70 

24 

B 

5.97 

40 

23. GO 

111.7 

3.G7 

7 

A 

5.99 

40 

33.79 

114.6 

5.07 

32 

B 

5.99 

40 

33.75 

115.5 

4.99 

9 

A 

C.00 

40 

43.04 

112.7 

G.G4 

15 

A 

5.99 

30 

33.79 

114. G 

5.07 

19 

A 

5.98 

20 

33.52 

114.1 

5.0G 

30 

B 

5.98 

20 

34.04 

115.4 

5.03 


Table 2 Test conditions for Rev— C configuration. 


Run 

Moc 

o 

(log 

Px 

lb/ft- 

Too 

(log R 

Re/ 

x 10"“ 

G 

5.9G 

30 

5.7G 

109.4 

0.93 

5 

5.97 

30 

11.80 

112.0 

1.83 

4 

5.99 

30 

23.35 

111.8 

3.G4 

7 

G.00 

30 

33.56 

114.5 

5.05 

31 

5.99 

20 

23.04 

111.1 

3.63 

2 

5.99 

40 

23.29 

111.4 

3.G5 


fluorescence intensity images of an illuminated model 
exposed to the wind-tunnel flow, surface temperature 
maps can be calculated on portions of the model in 
the camera’s view. A temperature calibration of the 
system conducted prior to the test provides data nec- 
essary to convert the two-color images to temperature. 
Acquiring images (temperatures) at different times in 
the wind-tunnel run enables local heat transfer coef- 
ficient to be computed. Comparison of heat transfer 
measurements obtained with the thermographic phos- 
phor technique to those obtained using conventional 
thin-filni resistance gauges have shown excellent agree- 
ment. 10 

In routine use of this technique, the phosphor 
coating (« 1 mil in thickness) has proven robust and 
does not typically require refurbishment, during a test. 
This technique, which has been widely used at the 
Langley Research Center, offers two distinct advan- 
tage's over conventional heat, transfer test methods. 
First , the measurements provide a quantitative resolu- 
tion of global surface' temperatures and [mating unlike 
discrete* gauge* measurements. In addition, the model 
construction, testing and data reduction can be per- 
former! much faster and cheapen- than other testing 
methods. 

Computational Methods 

Experimentally measured [mating data are* com- 
pared with predictions from a “benchmark" 1 Navier- 
Stokes code* (LAURA Langley Aerothernmdynamic 
Upwind Relaxation Algorithm), and an “engineer- 


ing” code* (LATCH Langley Approximate Three- 
Dimensional Convective Heating), which is based on 
a boundarv-layer-inviscid method. 

LAURA Code 

The LAURA code* 11 l;i is a three-dimensional, 
finite-volume. Navier-Stokes code* based on an upwind 
relaxation algorithm that can compute both perfect 
and real gas flows. Inviscid fluxes are computed us- 
ing Roe's averaging 1 ' 1 and Yee*s 1;> Symmetric Total 
Variation Diminishing (ST YD) to achieve second or- 
der accuracy away from discontinuities. Second order 
accurate, central differences are used to compute the 
viscous flux. 

The governing equations are solved by march- 
ing in ‘'pseudo' 1 time until a steady state solution is 
achieved. The treatment, of the governing equations 
is described as point implicit because variables at the 
cell center are treated implicitly, whereas the latest 
available data arc* used for the other terms. With this 
strategy, updating re 11 centered variables requires the 
inversion of a much smaller matrix and only one level 
of storage. The LAURA code has been shown to coin- 
pan* well with both flight 10 and wind tunnel 1 ' data. 

LATCH Code 

The LATCH code 18, 19 is an approximate three- 
dimensional heating cock* based on the axisymmetric 
analog for general three-dimensional boundary lay- 
ers." 0 In this method the boundary-layer equations are 
first written in a streamline-oriented coordinate sys- 
tem (a ,/Un) where a is the distance measured along 
an inviscid streamline. A is tangent to the surface and 
normal to the streamline direction, and n is normal 
to the surface. If the viscous crossflow in the bound- 
ary layer is small and can be neglected (as it can be 
when the streamline curvature is small or when the 
wall is cold), the boundary-layer equations reduce to 
an axisymmetric form if a is interpreted as the dis- 
tance* along an “equivalent" axisymmetric body and 
the metric coefficient, h, associated with the spread- 
ing of the streamlines, is interpreted as the radius 
of the equivalent axisymmetric body. This greatly 
simplifies the boundary-layer solution and means that 
approximate three-dimensional heating rates can be 
computed along an individual streamline independent, 
of other streamlines using any axisymmetric heating 
prediction method. The approach is simplified further 
by using an approximate integral heating technique,” 1 
which has been shown to agree with more detailed 
finite-difference boundary-layer solutions, to calculate 
heating rates. In practice, solutions are carried out 
along multiple streamlines simultaneously in a march- 
ing fashion so that the approximate boundary solution 
over the complete body is obtained in a single pass 
down the body requiring about 5 10 minutes on a 
desktop workstation . 

To apply this method, the inviscid streamlines 
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and metric coefficients are required and must he ob- 
tained from a three-dimensional inviscid solution. In 
this paper the inviscid solutions were obtained using 
the inviscid version of the LAUR A code, described pre- 
viously, and the inviscid version of the DPLUR code, 22 
a flowfield code that has been optimized for running 
on parallel computers. These solutions require more 
time than the LATCH boundary-layer solutions, but 
much less time than the full Navier-Stokes solutions. 
The primary advantage of the LATCH code is that 
reasonably accurate heating solutions can be obtained 
in much less time than for a Navier-Stokes code and 
thus more solutions can be computed for the same ex- 
penditure of resources. The relatively fast turnaround 
makes LATCH an extremely useful code for vehicle de- 
sign. The LATCH code has been shown to compare 
well with both wind tunnel and flight data and with 
other more detailed fiowfield codes. 19,23 

All of the calculations for wind tunnel conditions 
presented in this paper for both the LALIRA and 
LATCH codes were run for a constant wall temper- 
ature of 540 deg. R (80 deg. F). The heating data for 
these conditions are presented in terms of heat transfer 
coefficient to minimize the effect of wall temperature 
on the results. The flight calculations were run at the 
local radiation equilibrium wall temperature assuming 
a surface einmisivity of 0.9. 

Results and Discussion 

In this section, heating predictions from the 
LATCH and LAURA code are first compared with 
experimental heating measurements obtained in the 
20-Inch Mach 6 Tunnel on a D Loft forebody model 
and a Rev-C “tip-to-taiF model. Predictions from the 
LATCH and LAUR A codes are then compared for the 
peak laminar heating point on the X-33 design trajec- 
tory. 

D— Loft Forebody 

As a part of the X-33 boundary- layer transition 
study, 5 a 10-inch forebody model of the D-Loft config- 
uration (approximately 0.019 scale) was constructed. 
This model represents approximately the forward 70 
percent of the complete configuration from the nose 
to just ahead of the canted-fin, fuselage juncture (see 
Figs. 1 and 2). A perspective view of the forebody 
configuration is shown in Fig. 3. The model size rep- 
resented a compromise between the desire to test the 
largest possible model to enhance the probability of 
obtaining “smooth body or natural’' transition over a 
range of angles of attack (20° < a < 40°) and the 
desire to minimize the disturbance or blockage of the 
flow in the tunnel. 

Although no experimental hypersonic pressure 
data have been obtained on any of the X-33 config- 
urations, calculated results are available. The calcu- 
lated surface pressure distributions along the wind- 



Fig. 3 X-33 D— Loft forebody configuration 

ward symmetry plane are presented in Fig. 4 for 
a — 20°, 30°, and 40° angle of attack. Only the in- 
viscid results are presented since both the inviscid and 
viscous LAURA calculations are in good agreement. 
The pressures increase with angle of attack as would 
be expected. There are two regions of over expan- 
sion and recompression: one near x/L « 0.1, and the 
second near x/L « 0.45. The strength of both reconi- 
pression regions increases with .angle of attack. For 
this Mach number, the nondimensional pressure at the 
sonic point is approximately 0.49. thus at a = 20" and 
30° the sonic* point is on the nose ahead of the first re- 
compression but at a — 40° the sonic point is located 
much further downstream at x/L ~ 0.35. 



Fig. 4 Windward symmetry plane inviscid pres- 
sures on the D— Loft forebody configuration at 

Mo C = 6. 

Inviscid surface streamlines oil the lower surface 
are presented in Figs. 5 -7 for a = 20.° 30,° and 40.°, 
respectively. At a = 40,° the pressure 1 on the symme- 
try plane is higher than for the lower angles of attack 
and the streamline patterns indicate that there is out- 
flow on the lower surface. However, at a = 20,° the 
windward symmetry plane pressure is lower and the 
streainlne patterns indicate much more inflow than at 
the higher angles of attack. 
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Fig. 5 Inviscid surface streamlines on windward 
side the D-Loft forebody configuration at M ^ = 6, 

n = 20". 




Fig. 6 Inviscid surface streamlines on windward 
side the D-Loft forebody configuration at il/oc = 6, 
a = 30" . 




Fig. 7 Inviscid surface streamlines on windward 
side the D-Loft forebody configuration at = 6, 

n = 40". 


Windward symmetry plane heating results at o = 
40° and i?r/, - 1.83.rl() fi , an 1 presented in Fig. 8. 



Fig. 8 Windward symmetry heating on the D-Loft 
forebody at 4/^ = G, o = 40°, and Rcj. — 1.9^10*'. 

The heating results an* presented as a heat transfer 
coefficient ratio h/hy/i versus x/L. where h is the lo- 
cal heat transfer coefficient and hyy is the stagnation 
point heat transfer coefficient calculated using the the- 
ory of Fay and Riddell 24 for a sphere with the radius 
approximately equal to the nose radius of the model 
(R n = 1.033 inches). Presenting the heat transfer re- 
sults in this form (i. e. h/hyn) has two benefits. First, 
the sensitivity of the results to wall temperature vari- 
ations is essentially removed. This is important for 
these wind tunnel data because the wall to total en- 
thalpy ratio is relatively high (h u /h fl(r % 0.6) and the 
measured heating rate 1 is significantly affected by vari- 
ation in the wall temperature 1 while the heat transfer 
coefficient is not. For the 1 liypeTSonic portion of the 
X-33 flight trajectory the wall to total enthalpy ratio 
is lower (0.2 0.3) and the effect of wall temperature on 
heat. transfer rate is much lower than for the Mach 6 
wind tunned. The second benefit is that, for the range 1 
of Reynolds numbers in the 20-Inch Mach 6 tunnel, 
this ratio is relatively insensitive 1 to Reynolds number 
variation for laminar flow because 1 both the numerator 
and denominator have 1 a similar dependence. 

In figure 1 8. both the LATCH and LAURA cal- 
culations are 1 compared with the experimental heating 
data along the windward symmetry plane for an an- 
gle 1 of attack of 40" and a Reynolds number based 
oil model length of 1.9,rl0 () . In both sets of exhibi- 
tions the flow has been assumed to be 1 laminar because 
of the relatively low Reynolds number for this case. 
The LAURA predictions are in good agreement with 
the measured heating over the 1 entire 1 length of the 
model. The 1 differences are 1 largest in the 1 stagnation ro- 
gion whore the measured heating is approximately 10 
percent higher than the 1 LAURA prediction but down- 
stream of the stgnation region the LAI R A predictions 
generally agree with the measured heating data within 
about 5 percent or levss. The LATCH predictions are 1 in 
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good agreement with both the LAURA predictions and 
with the measured heating data over the entire length 
of the model. The good agreement between LATCH 
and the other results for this case is encouraging. 

For this same case, lateral distributions at four 
axial stations x/L — 0.3, 0.5. 0.7, and 0.9 are presented 
in Figs. 9 12. respectively. The heating results are 
plotted in the form of h/hfrR versus y/L where y is 
the lateral coordinate on the body which is zero in the 
symmetry plane and L is the model length (L = 10 
inches). 

At x/L = 0.3 (Fig. 9), the heating is a minimum 
near the symmetry plane (y/L = 0) and rises to a. peak 
in the chine region. Both LATCH and LAUR A pre- 
dict an increase in heating in the chine region as would 
be expected, with LAURA predictions being approx- 
imately 10 percent higher than LATCH predictions. 
However, both are significantly lower than the mea- 
sured data in this region (25-30 percent). This was 
not expected since the flow was thought to be lami- 
nar at this Reynolds number, and based on experience 
with other configurations such as the Shuttle Orbiter 
(which has a much smaller chine radii), it would be 
expected that the predictions would be in better agree- 
ment with laminar data. The LAURA solution was 
run for a sufficient time to achieve convergence for the 
grid that was used (1G1 x 181 x 64). A grid conver- 
gence study is planned but has not been completed. 
However, based on LAURA calculations on other con- 
figurations, such as the Shuttle Orbiter, it is thought 
that the grid is sufficient to produce accurate heating 
predictions. Similar results are noted for the other lat- 
eral heating distributions at x/L = 0.5, 0.7. and 0.9 in 
Figs. 10, 11, and 12. 



Fig. 9 Lateral heating distribution on the D— Loft 
forebody at x/L = 0.3, — 6, a = 40°, and Re l — 

1.9xl0 6 . 

For this angle of attack (a — 40°), heating distri- 
butions along the windward symmetry plane for three 
higher Reynolds numbers, Rei = 3.7, 5.0, and 6.6 
x 10, 6 are presented in Figs. 13, 14, and 15. For 
Re L = 3.7./T0 6 (Fig. 13), there is good qualitative 



Fig. 10 Lateral heating distribution on the D— 
Loft forebody at x/L = 0.5, Aloe =6, a = 40°, and 
Rc l = 1.9x10*. 



Fig. 11 Lateral heating distribution on the D— 
Loft forebody at x/L = 0.7, A/*, = G, o = 40°, and 
Re L = 1.9x10*. 



Fig. 12 Lateral heating distribution on the D- 
Loft forebody at x/L = 0.9, = 6, a = 40°, and 

Re L = 1.9x10*. 
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agreement. between the measured heating data and the 
laminar calculations of both LATCH and LAURA un- 
til near x/L % 0.6 where tin* measured heating rises 
rapidly indicating that boundary layer undergoes tran- 
sition. The heating continues to rise downstream until 
it exceeds the turbulent LATCH calculations near the 
end of the body suggesting that the How is close to 
fully turbulent, in this region. Experimental data from 
two different runs have been included in the figure. 
The two experimental heating results are in reason- 
ably good agreement.. 

For 7?c/ = 5.0*10° (Fig. 14). the measured heat- 
ing data follows the laminar LATCH calculations to 
a value of x/L % 0.1 where it begins to depart from 
the prediction. No LAURA calculations are available 
at this Reynolds number. At x/L « 0.5 -0.6, the mea- 
sured heating rises rapidly to a value slightly above the 
turbulent LATCH prediction and maintains that level 
until near the end of the body. In general the turbu- 
lent LATCH predictions agree with the measured data 
to within 10 12 percent. 

For the highest Reynolds number tested. R< 7, = 
G.G.rl0 (l (Fig. 15), laminar predictions are available 1 
from both LATCH and LAURA. The trend of the ex- 
perimental and predicted results are similar to those at 
Rx 1 = 5.0x10° (Fig. 14). except that the departure be- 
tween the measured heating and the* predictions near 
x/L % 0.1 is much more* apparent. It is interesting 
to note* that the location where the measured heating 
departs from the laminar predictions is just down- 
stream of the* first roeompression region observed in 
tin* predicted pressure distributions (Fig. 4). It may be 
possible' that this roeompression is introducing a dis- 
turbance 1 in the flow that causes the heating to depart 
from the expected laminar trend and eventually builds 
up and produce's the classical “transition” behavior 
that is observed further down the body. It may also be 
possible that there is sufficient, surface* roughness in the 
nose region to produce a disturbance in the* flow that 
causes the downstream heating to increase above the 
laminar level until it has grown sufficiently strong to 
promote the classical transition behavior observed fur- 
ther downstream. A behavior similar to this has been 
observed by Berry, et al. 25 during tests of discrete sur- 
face roughness on the Shuttle Orbiter. Whatever the 
cause, the heating aft of x/L ~ 0.1 does not follow the 
expected laminar trends. The heating when plotted 
as h/hh ff should remain approximately constant at 
a fixed location on tin* model with increasing Reynolds 
number for laminar flow.. Comparing results for the 
various Reynolds numbers. Figs. 8 and 13 15, it can be 
seen that this is the observed behavior for the laminar 
calculations but it is not for the experimental measure- 
ments: thus, there must be some kind of disturbance 
that is causing the flow to depart from laminar. If this 
is true it may also be affecting the lateral heating in 
the chine region that was observed previously at the 


lower Reynolds number in Figs. 9-12. 



Fig. 13 Windward symmetry heating on the D— 
Loft forebody at M ^ = G, <\ = 40°, and Rri = 
3.7*10“. 



Fig. 14 Windward symmetry heating on the D— 
Loft forebody at Mol = G, a = 40", and /?< /. = 
5.0x10“. 

For the highest. Reynolds number case. 7?c/ = 
6.6x10°, lateral distributions at the same* four axial 
stations, x/L — 0.3. 0.5. 0.7, and 0.9, as shown pre- 
viously at a lower Reynolds number are* presented in 
Figs. 16 19. From results presented previously, (see 
Fig. 15), it was shown that on the windward syinine- 
trv plane* the measured heating data departed from 
a laminar behavior downstream of x/L % 0.1. A 
similar behavior in the lateral heating distributions 
is observed in Figs. 16-19. The difference* between 
measured heating and laminar predictions increase's 
with each successive downstream lateral cut until at 
x/L = 0.7 (Fig. 18), the* heating on the center por- 
tion of the model agrees with the* turbulent LATCH 
predictions. The* spike in heating near y/L % -0.2 ob- 
served in this cut plane is caused by a small damaged 
area on the* model due to a particle in the flow strik- 
ing the phosphor coating. This unintended roughness 
dement caused localized early boundary transition. 
At x/L = 0.9, the localized transition region caused 
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Fig. 15 Windward symmetry heating on the D— 
Loft forebody at M oo = 6, o = 40°, and Ret = 6.6x1 0 6 . 

by the “trip” has merged with the wedge of turbu- 
lent flow over the central portion of the model and 
the measured data are in good agreement (10 percent) 
with the LATCH turbulent predictions. In fact, for all 
cases where the flow was clearly turbulent, the LATCH 
calculations agree well with the experimental measure- 
ments. 



Fig. 16 Lateral heating distribution on the D- 
Loft forebody at x/L = 0.3, Moo = 6, « = 40", and 
Rci = 6.6xlO fi . 

Finally, for a Ret « 5x1 0°, windward symmetry 
plane heating distributions for this configuration are 
presented for two lower angles of attack (a = 20° and 
30°) in Figs. 20 and 21. 

For a = 20° (Fig. 20), two sets of data obtained 
on different models (see Table 2), are included. Both 
sets of measurements are in excellent agreement with 
each other and also with laminar heating computed 
by LATCH. For this case the flow appears to be com- 
pletely laminar. 

At a — 30° (Fig. 21), the measured heating agrees 
with the LATCH laminar predictions for x/L < 0.6 
and then rises rapidly, indicative of transition, to a 
level slightly above the turbulent LATCH prediction. 
These results (shown in Figs. 20, 21, and 14) graphi- 

8 



Fig. 17 Lateral heating distribution on the D— 
Loft forebody at x/L = 0.5, M 0 c = 6, o = 40°, and 
R( l ~ 6.6xl0 6 . 



Fig. 18 Lateral heating distribution on the D- 
Loft forebody at x/L = 0.7, A/oe = 6, a = 40°, and 
Rci = 6.6xl0 6 . 



Fig. 19 Lateral heating distribution on the D- 
Loft forebody at x/L = 0.9, A/oc = 6, a = 40°, and 
Re i = 6.6xl0 6 . 
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Fig. 20 Windward symmetry heating on the D- 
Loft forebody at M x = G, n = 20°, and Ret — 
5.0;r 10°. 

cally illustrate tin* forward movement. transition with 
increasing angle of attack discussed by Thompson, 
et al.‘* 



Fig. 21 Windward symmetry heating on the D- 
Loft forebody at = G, a = 30°, and Rcl = 

5.0*10*. 

Rev-C Configuration 

Heat transfer tests have been conducted on a Rev 
C tip-to-tail model that had a length of 10 inches from 
the tip of the nose to the end of the engine block (ap- 
proximately 0.0132 scale). In this section measured 
heating data from these tests an 1 presented along with 
predictions from the LATCH code. 

First, heating distributions along the windward 
symmetry plane for Rcl ~ 3.G*10 6 and a = 20°, 
30° and 40° an 1 presented in Figs. 22. 23. and 24. 
respectively. In general, both the measured heating 
and the LATCH predictions increase with angle of 
attack would be expected. However, for all angles 
of attack, the experimental data are higher than the 
LATCH calculations in the nose region, x/L < 0.3, 
and somewhat lower for x/L > 0.3. The reason for 
this disagreement is unknown. Although, no LAURA 


calculations art* available for this configuration at these 
conditions, based on the comparisons between LAURA 
and LATCH for the D Loft forebody model shown 
previously and on other comparisons at flight condi- 
tions (which will be shown in the next section) LAURA 
and LATCH would be expected to give similar predic- 
tions for this case. However, as pointed out previously, 
the models have not been tested for quality assurance* 
and differences in shape between the wind tunnel and 
computational model geometries could cause some of 
the observed difference* between experiment, and pre- 
diction. For all future* X-33 heating tests, quality 
assurance measurements have been planned as an in- 
tegral part of the model construction process. 



Fig. 22 Windward symmetry plane heating on 
the Rev-C configuration at il/ x — G, <\ = 20°, and 
R<l = 3.6*10*. 



Fig. 23 Windward symmetry plane heating on 
the Rev— C configuration at i\/ x = G, o = 30", and 
Rcl = 3.6*10*. 

For n = 30° and Rcl = 3.6* 1() 6 , lateral heat- 
ing distributions at x/L = 0.3. 0.5, 0.7. 0.8. and 0.9 are 
presented in Figs. 25 29. respectively. The two most 
forward stations at x/L = 0.3 and 0.5 (Figs. 25 and 2G) 
are ahead of the* canted-fin. fuselage juncture* and the 
heating distribution is similar to that shown previously 
for the D Loft forebody mod(*l. Both the experimen- 
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Fig. 24 Windward symmetry plane heating on 
the Rev-C configuration at = 6, a = 40°, and 
Re L = 3.6*10°. 


tal data and the LATCH computations are lowest near 
the symmetry plane ( y/L = 0) and increase to a maxi- 
mum in the chine region. The LATCH predictions are 
in reasonably good quantitative agreement with the 
measured heating near the symmetry plane but are 
significantly lower than the measured data in the chine 
region. This is similar to the results shown previously 
for the D-Loft. Additional heating tests are planned 
for the Rev-F configuration, which is a later version of 
the X-33 vehicle than those included in this paper. As 
a part, of these tests, extensive supporting calculations 
are planned to further investigate the consistent trend 
observed here of underpredicting the “laminar” chine 
heating. 

Further aft on the model, the lateral cuts at 
x/L = 0.7, 0.8 and 0.9 (Figs. 27, 28, and 29) in- 
clude data on the windward side of the canted fin. 
For x/L — 0.7 (Fig. 27) and 0.8 (Fig. 28). the heat- 
ing increases in a direction away from the symmetry 
plane (y/L — 0), reaches a peak in the chine re- 
gion (y/L % 0.2 -0.3), decreases around the side of the 
fuselage, and then increases again on the windward 
side of the canted fin. The LATCH predictions are 
in good qualitative agreement with the experimental 
measurements but again underpredict the heating on 
the chines. The overall comparison on the windward 
side of the canted fin is better than on the chines. For 
x/L = 0.9 (Fig. 29), the trends are similar to those 
at the two previous axial stations ( x/L = 0.7 and 0.8) 
hut lower heating is observed over the center portion 
of the model. This lower heating is caused by an ex- 
pansion region on the rearward portion of the lower 
surface. The tip of one of the canted was broken off 
during the expermental tests. This accounts for the 
asymmetry in the measured heating outboard of the 
chine region for negative values of y/L at this axial 
station (x/L = 0.9). 



Fig. 25 Lateral heating distribution on the Rev— 
C configuration at x/L = 0.3, M 0 c = 6, a = 30°, and 
Rc l = 3.6*10°. 



Fig. 26 Lateral heating distribution on the Rev- 
C configuration at x/L = 0.5, Af*, = 6, a = 30°, and 
Rc l = 3.6*10°. 



Fig. 27 Lateral heating distribution on the Rev— 
C configuration at x/L — 0.7, M 0 c = 6, o = 30°, and 
Re L = 3.6*10°. 
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Fig. 28 Lateral heating distribution on the Rev 
C configuration at x/L = 0.8, M^c = G, a = 30°, and 
Rvi. = 3.6x10“ . 


Lateral Heating on X-33 Rev-C Configuration 
x/L = 0.9 



y/L 


Fig. 29 Lateral heating distribution on the Rev 
C configuration at x/L = 0.9, Jl/oc = G, a = 30", and 
Rc r = 3.6x10“. 

Flight Predictions 

In this section computed heating results from the 
LATCH code are compared with computed heating re- 
sults from the LAURA code at the peak laminar heat- 
ing point on the nominal TPS (Thermal Protection 
System) design trajectory. The freest ream conditions 
for this case are given in Table' 3. 

Table 3 Freest ream conditions for flight calcula- 
tions 


Altitude = 

175.8 

kft 

M-x 

11.5 


a - 

3G.2 

deg 

Px = 

1.0636 

lb/ft.2 

lx 

47C.4 

deg R 

Re,. 

2.9 x 10“ 



Figure 30 shows a plan-form view of the lower 
surface of the X-33, D Loft configuration with com- 
parisons of the heating rate contours from the LATCH 
and LAURA cock's. The 1 LATCH results are shown on 


the bottom half of the figure and those from LAURA 
on the top. These' calculations were' for laminar flow 
and the wall temperature was assumed to be the* ra- 
diation equilibrium wall temperature 1 . The heating for 
this e*ase* is presented directly in tennis of heat trans- 
fen* rate (q, r ) because the wall to total enthalpy ratio 
is relatively low (0.2 0.3) and the heating rate is rela- 
tively insensitive to wall temperature. From Fig. 30, it 
is observed that the two codes are in good qualitative 
agreement for this case. The* heating rate contours in 
both the chine and cantod-fin regions are* very similar. 
Although the heating rate* contours near the windward 
symmetry plane* appear at first glance to be quite dif- 
ferent, the heating in that region is nearly constant, and 
the* two calculations are* actually in very good agree- 
ment. 



Fig. 30 Comparison of LATCH and LAURA heat- 
ing rate predictions on the lower surface of X-33 in 
flight. 

Heating distributions along the windward sym- 
metry plane* are* presented in Fig. 31 and as well as 
lateral distributions at several axial stations, r = 250, 
500. 550. 600. 650, and 700 iimhes. in Figs. 32 37. re- 
spectively. The* lateral cut. at x = 250 inches (Fig. 32) 
is on the forebody ahead of the canted fin. All of the* 
remaining lateral cuts (Figs. 33 37) are downstream 
of the the cant.ed-fin fuselage juncture* (see Figs. 1 and 
2). In general the LATCH and LAURA results are 
in very good agreement. (d= 10 percent) ove*r most of 
the windward surface* e*xcept in the chine region ne*ar 
the rear of the vehicle* where* the LATCH prcdictiems 
are lower than LAURA by approximately 20 30 per- 
cent.. The* agreement between LATCH and LAURA 
at this flight condition is encouraging because it sug- 
gests that the* LATCH code*, which has a much faster 
turn around time than LAURA, can be used to predict, 
flight heating rates for X-33 TPS design. 

Concluding Remarks 

Laminar and turbulent he*ating-rate* calculations 
from an “engineering" code (LATCH) and laminar 
calculations from a “benchmark’’ Navier-Stokes code 
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Fig. 31 Windward symmetry plane heating on the 
D-Loft configuration at Moc = 11.5, o = 36.2°, and 
Re- 1. = 2.9a:10 6 . 



Fig. 34 Lateral heating distribution on the D- 
Loft configuration at x = 550 inches, = 11.5, 

a = 36.2", and Rc,. = 2.9a-10 6 . 



Fig. 32 Lateral heating distribution on the D- 
Loft configuration at x = 250 inches, M 0 0 = 11.5, 
a = 3G.2°, and Re L = 2.9;rl0 6 . 



Fig. 35 Lateral heating distribution on the D- 
Loft configuration at i = 600 inches, Moo = 11.5, 
a = 36.2°, and Re L = 2.9xl0 6 . 



Fig. 33 Lateral heating distribution on the D Fig. 36 Lateral heating distribution on the D- 

Loft configuration at x = 500 inches, = 11.5, Loft configuration at x = 650 inches, Moc< = 11.5, 

a = 36.2°, and Re L = 2.9x10®. « = 36.2°, and Re L = 2.9.rl0 6 . 
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Fig. 37 Lateral heating distribution on the D- 
Loft configuration at r — 700 inches, M ^ = 11.5, 
n = 36.2°, and He /. = 2.9;rl0 <> . 

(LAURA) wen' ronipared with experimental wind- 
tunnel data obtained on several candidate configura- 
tions for the X-33 Phase II flight vehicle. The experi- 
mental data were obtained at a Mach number of G and 
a freestream Reynolds number ranging from approxi- 
mately about 1 x 10 (i to 8 x 10 (i /ft* Comparisons are 
presented along the windward symmetry plane and in 
a circumferential direction around the body at several 
axial stations at angles of attack from 20" to 40°. The 
experimental results include both laminar and turbu- 
lent how. 

For (\ — 40." the measured heating data along the 
windward symmetry plane exhibited a “noil-laminar 
behavior begining near x/L — 0.1 which caused the 
heating to increase above the laminar level long before 
‘■classical*' transition to turbulent, flow was observed 
and the departure from the laminar heating level in- 
creased with increasing Reynolds number. This be- 
havior was not observed at the lower angles of attack 
(i. e. n = 20" and 30"). The predicted inviscid pres- 
sure distribution along the windward symmetry plane 
exhibited two over expansion and recompression re- 
gions, one near the nose (x/L % 0.1) and the* other 
near the midpoint of the vehicle' (./*/ L « 0.45). These 
recompressions were observed to increase in strength 
with angle of attack. Disturbances emulating from the 
recompression in the nose 1 region may have disturbed 
the flow and contributed to the noil-laminar behavior 
of the* heating observed at the highest angle of attack. 

When the* flow was completely laminar, both 
code's predicted the measured heating along the 
windward symmetry plane* reasonably well, generally 
within 10 percent: however, both cock's under pre- 
dicted the' measured heating in the chine region, gen- 
erally by 25 30 percent.. When the flow was turbulent 
the LATCH code accurately predicted the measured 
heating rate's within about 10-12 percent. Both codes 
were used to calculate' heating rates over an X-33 vehi- 
cle at the peak heating point on the* design trajectory 


and they were* found to be in very good agreement 
except near the rear of the vehicle in the chine region. 
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